Direct Urea neutrino rate in colour superconducting quark matter 
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If deconfined quark matter exists inside compact stars, the primary cooling mechanism is neutrino 
radiation via the direct Urea processes d — > u + e + i? e and u + e — » d + v e . Below a critical 
temperature, T c , quark matter forms a colour superconductor, one possible manifestation of which is 
a condensate of (ud) quark Cooper pairs in an electric-charge neutralising background of electrons. 
We compute the neutrino emission rate from such a phase, including charged pair-breaking and 
recombination effects, and find that on a material temperature domain below T c the pairing-induced 
suppression of the neutrino emission rate is not uniformly exponential. If gapless modes are present 
in the condensed phase, the emissivity at low temperatures is moderately enhanced above that of 
completely unpaired matter. The importance of charged current pair-breaking processes for neutrino 
emission both in the fully gapped and partially gapped regimes is emphasised. 

PACS numbers: 97.60.Jd, 24.85. +p, 26.60.+C, 95.30.Cq 



The direct Urea neutrino emission process (in-medium 
/3-decay), which was originally introduced in studies of 
stellar collapse and supernova explosions [ll, plays a piv- 
otal role in the cooling of compact stars 0, 0j- When 
kinematically allowed it is one of the most efficient mech- 
anisms by which massive neutron stars composed only of 
ordinary nuclear matter can cool, allowing them to reach 
temperatures below the observational threshold within 
roughly a decade Q. 

It is possible that the proton concentration in the in- 
terior of such a star may be insufficient to satisfy energy- 
momentum conservation in the decay n — > p + e + D e . 
However, if the interior contains unpaired quark matter, 
then the reactions d — > u + e + v e and u + e — * d + v e can 
proceed easily since the quarks are ultrarclativistic, and 
interactions and/or finite quark masses guarantee that 
energy-momentum conservation is satisfied. The reac- 
tion rate is proportional to the strong coupling constant 
a s in the weak-coupling limit Q. These reactions pro- 
vide for prodigious neutrino emission, which dominates 
over all other sources of neutrino production. Hence, the 
presence of unpaired quark matter in a neutron star is 
expected to dramatically alter the star's cooling rate and 
consequently its temperature-age profile, in which case 
astronomical observations of this profile may be used to 
place constraints on the abundance of quark matter in 
stellar interiors. 

This straightforward connection may be undermined 
if quarks condense to form Cooper pairs in the dense 
stellar interior. (Aspects and consequences of such con- 
densation are described in Rcfs. In this instance, 
one might anticipate that the direct Urea neutrino rate 
at temperatures T < T c will be greatly suppressed owing 
to the cost of bridging the pairing gap A(T) via ther- 
mal excitation. (T c < 0.1 GeV is the temperature above 
which the pairing gap disappears.) 

Herein we show that this is not necessarily the case be- 
cause at moderate temperatures, T < T c , modifications 
of the phase-space occupation factors resulting from the 



gap are non-exponential, and the breaking and forma- 
tion of pairs also contributes to the neutrino rate. In the 
following we describe a calculation that illustrates these 
features and supports a view that the suppression just be- 
low T c is not uniformly exponential. This enables us to 
quantify the direct Urea neutrino rate in superconduct- 
ing quark matter, which may be relevant to the cooling 
of compact stars. 

The central densities of neutron stars may be as much 
as (5 — 10)po depending on the underlying equation of 
state of dense matter; here po = 0.16 fm -3 is the nuclear 
saturation density . At such densities an ultrarelativis- 
tic two-flavour quark gas is characterised by a chemical 
potential p U: d ~ 0.4 — 0.5 GeV > Aqcd- Perturbation 
theory is not a valid tool on this domain of \x u ^d and 
hence the true nature of the ground state of quark mat- 
ter is thereupon uncertain. Far from this scale, namely, 
at [i u _d 3> Aqcd, perturbation theory predicts that quark 
matter exists in a colour-flavour-locked phase (CFL) 0. 
However, that is only one amongst many phases which 
may be realised as the chemical potential is reduced from 
asymptotic values to that characteristic of a neutron star 
interior (NB. Measured in scales typical of QCD, 
neutron stars are cool; viz., T ~ 0.) 

We will explore a minimal realistic scenario of two- 
flavour (u, d) quark matter, with spin-singlet chan- 
nel pairing, constrained by charge neutrality and (3- 
cquilibrium. Since increasing \i promotes quark-quark 
pairing, while increasing quark mass opposes it; viz., the 
magnitude of any putative gap, A(T), increases with in- 
creasing chemical potential but decreases with increas- 
ing mass, this picture can be valid when the relevant 
s-quark mass-scale @ satisfies MKpl/T) > ^A(T). For 
p ~ 0.4 GeV and A(T) < 0.1 GeV [ljj this means our 
scenario may be valid for M s (p,,T) > 0.2 GeV. On this 
domain the s-quark is too heavy to influence pairing be- 
tween the light-quarks. Hence it plays no role in our 
subsequent analysis. We note that M s (fi ~ 0.4 GeV, T < 
T c ) > 0.2 GeV is plausible |n|. 
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It is plain that the enforcement of charge neutrality and 
/3-cquilibrium entails a mismatch between the u- and d- 
quark chemical potentials; i.e., Sfi ~ (/i^ — Mu)/2 = /i e /2. 
The neutrino radiation rate therefore depends on the ra- 
tio £ := A/Sfi. On the domain £ > 1 the relevant 
quasiparticle dispersion law is gapped for all momen- 
tum modes. However, for £ < 1 there is a measurable 
domain of ungapped modes |l2j . Neutrino radiation in 
this regime has previously been studied for three flavour 
(u, d, s) pairing |l3|. NB. A nonzero /3-decay matrix ele- 
ment is only possible with a dressed-quark dispersion law. 
Such modifications can arise through, e.g.: pairing; inter- 
action effects that are linear in a s ; a dynamically gener- 
ated quark mass 0] ; and non- Fermi-liquid behaviour 01 • 
In the following we consider massless quarks, interaction 
corrections that are linear in a s , and the influence of 
pairing. 

Cooper pair breaking and formation processes play a 
central role during the neutrino radiation era both for 
superfluid neutron stars 0] and compact stars featuring 
colour superconducting quark cores 13, 161. The inclu- 
sion of pair-breaking and recombination effects below T c 
is therefore an important feature of our study. The micro- 
scopic origin of the Cooper pair breaking and formation 
processes and their contribution to the neutrino emission 
from neutron star matter was explored in a number of 
works 0, 0, 0, 1^3 ■ The counterparts of these pro- 
cesses in quark matter, i.e. the case where the neutrino 
emission is driven by weak neutral current interaction 
between quarks and neutrinos, was studied in Ref. |2ll| . 
Herein we explore the role of the weak charge current 
pair-breaking process in neutrino radiation, which was 
found to be important in the nuclear Urea process |l9j | . 

In order to calculate the total neutrino emission rate 
we note that, since the rates for /3-decay and electron 
capture are identical in /3-equilibrated matter, it is only 
necessary to compute one rate, say d — > u + e + z? e , and 
subsequently multiply the result by two. We employ the 
Kadanoff-Baym formalism for neutrino transport [23.l23| . 
wherein the neutrino emissivity, defined as the neutrino 
rate per unit volume, can be related to the IF-boson 
polarisation tensor in the superconducting phase; namely 
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FIG. 1: One-loop V^-polarisation tensor: (a) normal piece, 
which is the sole contribution above T c ; and (b) "anomalous" 
piece, which is proportional to A 2 (T) and vanishes at T c . The 
imaginary part of (b) (nonzero below T c ) contributes the v- 
rate from charged-pair breaking and formation. 



not have support for £q < 0). The lepton tensor in Eq. (JIJ 
is 

A^( qi ,q 2 ) = Tr [y (1 - 7 5 ) rf l7 *(l - 7 5 ) 4 2 \ 
= 8^2 + ?2% A -91-92 g» X + i^ aX(i qi a q2 P ] ■ (2) 

At one-loop level, the retarded M^-boson polarisation 
tensor in the superconducting phase, n^ A (fc), is obtained 
by combining the two physically distinct contributions il- 
lustrated in Fig.^ To evaluate this tensor, we require ex- 
pressions for the quark propagators in the two-flavour su- 
perconductor. These are detailed in [24|, where the gluon 
polarisation was studied. For our purpose, we replace the 
gluon vertex by the weak vertex r± = 7 M (1 — 75) ® t±, 
where t± = (t\ ± t 2 )/2 are flavour-raising and lowering 
operators constructed from Pauli matrices. In addition, 
antiparticle contributions are neglected because we work 
at large quark chemical potential. 

It follows that diagram (a) in Fig. ^ represents 
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with 



H 



(«) 



Tr[ 7/ ,(1 - 75)A+(p) 7 o7a(1 - 75)A+(k) 7o ] , (4) 



where the trace is over colour, flavour and spinor indices 
and A ± (k) = (1 ± 7 • k)/2. In Eq. J3J, = k - ft 

and E~ = p - p., with p = (p d + p, u )/2, = 



x n F {-uj 1 )n F (-u 2 )K^ x (q u q 2 ) SJmn^ A (^) , (1) \J E^J + A 2 (T), and the Matsubara sum is over dis- 



wherc: G — Gp cos 9c, with Gp being the Fermi weak 
coupling constant and Be the Cabibbo angle (cosflp = 
0.973); (71,2 = (oJi t2 ,q 1 2 ) arc the 4-momenta of the lep- 
tons (e and v, both considered massless and labelled, 
respectively, by indices 1 and 2) and I is the VF-boson 
momentum; J := J d 3 qi/[(2ir) 3 2u>i] are phase-space in- 
tegrals; is the Bose distribution function for the 
IF-boson energy; and np(— Wi, 2 ) are the Pauli blocking 
factors for electrons and neutrinos in the final state. (NB. 
Emission is characterised by loi_ 2 > and hence the Dirac 
delta function in Eq. ensures that the integral does 



crete frequencies po n = (2n + l)inT, with ko n = Pon +£o- 
We observe that in the normal phase, defined by A = 0, 

A^ X H^ =64£JiW 2 (l-qi-k)(l-q 2 -p) =0, unless the 
quark masses are nonzero or one includes perturbative 
corrections to the free-quark dispersion relation. This fol- 
lows from four-momentum conservation, which for mass- 
less free-quarks requires all four particles to be collincar; 
viz., qi • k = 1 = q 2 • p. Once interactions modify the 
quark dispersion relation, this is no longer true, and Q 



A^> = 64^ 2 ^ 



(5) 
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provides the leading-order pcrturbative contribution 
from diagram (a) to e„ in Eq. |T|I, 

Turning to the anomalous contribution, (b) in Fig. Q 
which has gap insertions in the quark lines, we obtain 



1 



5>, 



0) 

A [(k 0n + 5^ - # 2 ][(po„ + SnY - if] ' 



(6) 



with 



#2 = Tr Ml ~ 75)A + (p) 7 57a(1 + 7 5 )A + (k) 75 ], (7) 



where = \J + A 2 (T), with E+ = k + p.. In 
Eqs. Q and ©, the trace over the colour-flavour struc- 
ture of the condensate yields a factor of two. 

The material features of the anomalous contribution, 
Eq. 10, are that it appears only in the superconduct- 
ing state and vanishes in the limit T — > T~ . For our 
argument and illustration it is sufficient to employ a pa- 
rameterization of the gan 1251 that is based on the BCS 
model calculation of Ref . Ufa : 



A(T) = A(0) 1 - (T/T c ) 



\3.4 



0.53 



0(T C - T) . (8) 



While modifications to this profile have a modest quanti- 
tative impact on our results, they are immaterial to our 
primary conclusions. 

Since we have simple algebraic expressions for each el- 
ement that appears in Eqs. (jSJ and J^J, it is straightfor- 
ward to evaluate the Matsubara sums and subsequently 
to obtain the imaginary part of the retarded polarisation 
tensor via analytic continuation |27j| : i.e., £o ~ * ^o + *£■ 
From Eq. © one obtains the normal contribution: 



Smn fl ; (a) (i) = -- 



d3 P rr{a) 

(2tt) 3 ^ a 



x [(f P ' + - w - ) i^+i + 1) (£p + 1) 

- (i - - / P ~+7) W) + 1) (£p - 1) 

- (fp- /;+!) <w + ) ( £ P+e - i)> (£, - 1)] > (9) 

where f^' a = + aSp,), r = p,p + £, are equilib- 

rium free particle Fermi distribution functions, <5„(£q^) = 
6(£o + a£- +K p+l + 25^) and £f = E±/£±. The anoma- 
lous contribution is obtained from Eq. © 



p+>- 



(2tt) 3 



ma t- t+ 



x [(/,"•" - ~) - (/p+7 - fp- + ) W + ) 

- (i - - W)]> (io) 

where «y o (^) = 5(4 + a£ + /^+ +£ ). 

The Dirac-delta functions in Eqs. J^J and (|1(J|> make 
the kinematics of scattering and pair breaking transpar- 
ent. The 5(Iq~) and S{Iq + ) terms describe scattering of 



u- and d-quark quasiparticles from the condensate. The 
S(Iq~) terms, on the other hand, express breaking of 
charged Cooper pairs when the total energy £q exceeds 
the pair threshold; viz., £$ > 2A. The term proportional 
to S(Iq + ) has been dropped as it does not contribute to 
the emission rate. The delta functions also simplify the 
evaluation of the integrals in Eqs. © and I|1U[) . More im- 
portantly, they place bounds on the magnitude of p 2 so 
that the integrals are finite. 

We now return to the neutrino emission rate of Eq. Q . 
In addition to Eq. (JSJ and at the same order we have 



A" A ffJ } = 32^ 2 



(1 + p • k) 



16a s 

~3tT 



pq2 



qiq 2 (l - p • k) - p • qi (1 + k • q 2 ; 



(11) 



Thus 



d i £ 

—— -S 4 (£- qi -q 2 )uj 2 {q 2 ) 

(2?T 4 ) 



with 



x n B (£ )n F (-uji)n F (-aj 2 )J : '{£o 1 £ 1 qi,q 2 ) , (12) 

T(£ Q ,£, qi , q 2 ) = A^ A ( qi , g 2 ) 3m n« A (£ Q ,£) 
= A^(q u q 2 ) [3mnJ W (f) + 9mny»(4 (13) 



It is straightforward to evaluate the neutrino emission 
rate in Eq. I|12|) . One of the integrations, say that over the 
neutrino's three-momentum, is made trivial by the factor 
S 3 (£ — qi — q 2 ): q 2 = i — q 1 , which leaves S(£o — cJi — uj 2 ) 
with u>2 = 1 92 1 = \£ — <?i|. This last distribution bounds 
the accessible domain of [q^; namely, £o — \£\ < 2|qr 1 1 < 
£q + \£\, and fixes the value of £ ■ q 1 within this domain. 
It remains to evaluate the integrals over \qi\, £ and £o, 
tasks we perform numerically, in the order listed. NB. 
The Bose factor n B (£ ) guarantees rapid convergence of 
the £o integral. 

Our calculated result for the emissivity of superfluid 
quark matter at T = T c matches that of the unpaired 
state 0: 



e = (914/315) G F a s fi d [i u fi e T 6 



(14) 



to within a few percent. The small deviation owes to our 
more complete treatment of nonzero temperature effects. 

In Fig.0 we depict the ratio e^/eo as a function of 
t = T/T c . The upper panel makes plain that neutrino 
production is suppressed in the presence of complete pair- 
ing; i.e., when all momentum modes are gapped. How- 
ever, when the magnitude of the gap is moderate, that 
suppression is not uniformly exponential. Rather, it is 
approximately linear on a material domain. For the pa- 
rameters we have chosen, that domain is 0.4 < T/T c < 1. 
Naturally, at any given value of t the emissivity decreases 
with increasing Q: a larger gap must be overcome. In the 
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FIG. 2: Temperature dependence of the emissivity Eq. (|12|l 
normalized to its value at the critical temperature. Upper 
panel - fully gapped regime, £ > 1; and lower panel - partially 
gapped regime, ( < 1. In both panels the dashed and dash- 
dotted curves are, respectively, the normal and anomalous 
contributions; and the solid line is their sum. In the lower 
panel, for reference, the horizontal line e v /ta — 1 is the result 
for unpaired matter. The curves were obtained with /j,d = 400, 
flu = 300 and /j, e — Hd — fi u = 100 MeV, values which are 
typical of two-flavour colour superconducting quark matter. 

fully-gapped regime, for straightforward kinematic rea- 
sons, the anomalous contribution vanishes in the limits 
t -> and t -> 1. 

The emissivity in the partially-gapped regime is il- 
lustrated in the lower panel of Fig.0 V£ € (0,1) the 
neutrino emission is enhanced cf. Eq. I|14fl, The normal 
contribution is magnified because the gap persists for 
a measurable set of momentum modes and hence more 
phase space is accessible to these modes than would oth- 
erwise be the case. The effect disappears as £ — > 0. The 
behaviour of the anomalous contribution in this regime 
is also understandable: it is nonzero but decreases uni- 
formly in magnitude as C —> 0, driven by the reduction 
in phase space and the overall multiplicative factor of A 2 
in Eq. 10 . The result for unpaired matter is recovered at 
t = 0. 

In the neighbourhood of t = 

e„/eo oc {(7(C) + exp[(C - l)<WT]}-\ (15) 

where C(£) ~ 0(1) is independent of temperature; viz., 
in the fully-gapped regime the normalised emissivity is 
exponentially suppressed, while in the partially-gapped 
regime it tends to a constant. 



We have considered isospin singlet pairing and the be- 
haviour of e„ exhibited herein differs from that found in 
earlier studies of isospin-triplet pairing; viz., Refs. [2flL 
l3fj| |. which focused on reaction phase space, omitted the 
coherence factors and pair-breaking processes, and used 
a non-BCS dispersion relation; and Ref. [l9j , which in- 
cluded those effects along with a BCS dispersion relation. 

Our calculation in the two-flavour phase is indicative 
of similar behaviour in the neutrino rate below T c for 
the three- flavour superconductor (CFL phase). However, 
in that case quark quasiparticle excitations play a role 
secondary to CFL Goldstone bosons so far as neutrino 
emission and cooling of stars with a CFL core are con- 
cerned [jul, |33|. As a follow up to the present study, 
it may be useful to determine the neutrino opacity in a 
two-component quark superfluid. This, determined by 
scattering and absorption of neutrinos, might be accom- 
plished by following the pattern of Ref. [33. Neutrino 
emission via the Urea process from a spin-one colour su- 
perconductor has also been assessed |34l l35l |. Neutral 
current mediated pair-breaking processes would be rele- 
vant to spin-one phases of paired quark matter. 

We focused on charged-pair breaking because the neu- 
trino emission rate generated by direct Urea is far greater 
than that provided by the neutral analogue (neutrino 
bremsstrahlung). In addition, for the two- flavour case 
there are no charged Goldstone modes that can couple 
to neutrinos. The contribution to neutral-current neu- 
trino scattering and emission from such soft (lj <C A) 
collective excitations was highlighted for the CFL phase 
in Ref. [3j| . For example, their inclusion is essential to an 
internally consistent truncation of the Dyson-Schwingcr 
equations and the concomitant preservation of relevant 
Ward-Takahashi identities. As noted, such contributions 
are absent in the discussion of a two-flavour supercon- 
ductor and one can unambiguously obtain the transverse 
part of the W^-boson polarisation tensor without model- 
dependent assumptions about the dressed weak vertex. 

To close, we note that the temperatures within a star 
during the protoneutron phase of its evolution may sup- 
port a colour superconducting state in its dense core be- 
cause the critical temperature is anticipated to be large: 
T c x O.6A0 ~ 50 MeV. Neutrinos would interact often in- 
side quark matter and therefore serve as an electroweak 
probe of the superconducting phase. Numerical simula- 
tions of the early stages in the evolution of newly-born 
hot neutron stars indicate that these systems support a 
lepton fraction of « 0.4 (3(|. With the neutrinos trapped, 
this fraction represents electrons and neutrinos. 

Now, should a transition to neutral quark matter take 
place at this early stage, the system can satisfy the con- 
dition of /3-equilibrium (/x^ — fj, u = /j e — ) while sup- 
porting an electron fraction x e » 0.2. (NB. Neutrinos 
free stream in a cool star, in which case [l V( , ks and 
x e rs 0.01.) Furthermore, recent model estimates of the 
equation of state for /3-equilibrated charge-neutral matter 
suggest that in the dense core of a neutron star, despite 
its positive charge, two- or three-flavour Cooper-paired 
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quark matter can occupy a large fraction of the volume 
as one component of a mixed phase that also contains 
normal quark matter ■ 

Considered together with neutrino opacities, which are 
large at T ~ T c , it is apparent that the results illustrated 
by Fig.Elare pertinent to any detailed computation of the 
cooling history of a star that undergoes a phase transition 
from nuclear matter to superconducting quark matter. 
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